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Bubble Chambers provided the dominant particle detection technology in ac- 
celerator experiments for several decades, eventually falling into disuse with 
the advent of other techniques. We report here on the first period of op- 
eration of an ultra-clean, room-temperature bubble chamber containing 1.5 
kg of superheated CF3I, a target maximally sensitive to spin-dependent and 
-independent Weakly Interacting Massive Particle (WIMP) couplings. An ex- 
posure in excess of 250 kg-days is obtained, with a live-time fraction reaching 
80%. This illustrates the ability to employ bubble chambers in a new realm, 
the search for dark matter particles. Improved limits on the spin-dependent 
WIMP-proton scattering cross section are extracted from this first period. An 
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extreme intrinsic insensitivity to the backgrounds commonly limiting these ex- 
periments (a rejection factor for photon-induced electrons of ~ 10^°) has been 
measured in operating conditions leading to the detection of low-energy nu- 
clear recoils such as those expected from WIMPs. 

Introduction 

With approximately 85% of the total matter of the universe in a form which still eludes di- 
rect detection, the need for large-mass, background-insensitive detectors able to explore the 
very small couplings expected from Weakly Interacting Massive Particles (WIMPs) is urgent. 
While a member of this family of hypothetical particles, the Lightest Supersymmetric Partner 
(LSP, a.k.a. neutralino) stands out as one of the most likely candidates for the Dark Matter in 
galaxies (17]|2|), its predicted interaction rate via low-energy elastic scattering off nuclei is <<1 
event/kg-day for any target. In some unfavorable, yet entirely plausible models this can be- 
come a dismayingly small < 1 event/ton-year. To further aggravate this situation, the couplings 
still allowed by present direct searches are already so weak that the next generation of massive 
(O(lOO)kg) WIMP detectors will necessarily suffer from a penury of statistics in any dark mat- 
ter signal they may observe. For instance, unique signatures such as a ~5% annual modulation 
in WIMP interaction rate expected from the orbital motion of the Earth Q will soon require 
extreme exposures of O(100)Ton-year to become evident. Information about WIMP properties 
(e.g., mass of the particle) will be scarce. This clearly indicates the need for a variety of detec- 
tion techniques and targets, and for an emphasis within each experiment not only on improving 
sensitivity, but on developing schemes that allow the unequivocal identification of these parti- 
cles 0. 

Superheated liquids, in the form of bubble chambers, were extensively used as the detec- 
tion medium of choice in accelerator experiments spanning more than three decades. These 
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chambers led not just to numerous discoveries, but also to an effect nearly as important: an un- 
precedented ability to visualize fundamental particles and their interactions (|5]0). The advent of 
newer technologies led to their progressive obsolescence during the seventies. In a bubble cham- 
ber, the heat deposited by ionizing radiations along their path produces local nucleations of the 
vapor phase in a delicate (metastable) superheated liquid. Rapidly growing bubbles form along 
this path, are photographed, and the chamber is then reset by fast recompression to the stable liq- 
uid phase. Decompression follows, bringing the target below its vapor pressure at the operating 
temperature, sensitizing it to radiation and starting the cycle anew. In accelerator experiments 
it was sufficient to maintain the superheated state during the few milliseconds corresponding to 
the bunched arrival of incident particles. Uncontrollable boiling on porous materials (gaskets, 
metallic surfaces) prevented serious consideration of this type of detector in searches for rare 
events, where the time of arrival of the signal is not known. 

Taking several precautions towards the deactivation of inhomogeneous bubble nucleation 
centers, we have recently shown that it is possible to achieve what is in principle an indefi- 
nitely long stability in moderately superheated bubble chambers (jTD, thereby allowing their use 
in direct WIMP searches. In this realm of application, several advantages for this type of de- 
tector can be mentioned. Most important is the fact that the superheated liquid can be tuned 
to respond exclusively to particles having a large stopping power (dE / dx, energy loss per unit 
path distance). In this way muons, gamma rays, x-rays, beta particles, etc., all fall well below a 
bubble nucleation threshold, which is typically dE/dx > 50 keV/ fim during a WIMP search 
(Fig. 1). The target liquids are nevertheless sensitive to few-keV nuclear recoils such as those 
expected from the scattering of WIMPs, given their much denser energy deposition. The pro- 
cess of radiation-induced bubble nucleation is described within the framework of the classical 
"Hot Spike" model (8): for the phase transition to occur, the energy deposited by the particle 
must be larger than the thermodynamical work of formation of a critically-sized protobubble. 
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but this energy must also be lost over a distance comparable to the size of this protobubble, i.e., 
a minimum stopping power condition must be fulfilled. This condition leads to an advantageous 
insensitivity to the listed minimum ionizing backgrounds that normally plague WIMP searches 
(Figs. 1,2). Formal details on the theory of radiation-induced bubble formation in superheated 
fluids can be found in (|9]|73[77]) and references therein. As for the probability of spontaneous 
(homogeneous) bubble nucleation in the bulk of a superheated fluid, it only becomes sizeable 
a few degrees below the critical temperature of the target compound (|72)) . In normal operating 
conditions for a dark matter search this source of instability is entirely negligible. 

Experimental 

Chamber operation Early in 2005 we installed a chamber containing 1.5 kg of CF3I, an 
industrial refrigerant commonly used as a fire extinguisher, at the 350-foot depth of the NuMI 
(Neutrinos at the Main Injector) tunnel in Fermilab. Environmental neutrons can produce single 
bubbles such as those expected from WIMPs if they scatter just once in the target (Fig. 2). At this 
shallow depth, a combination of active and passive shielding can lead to a neutron background 
rate as low as 0.01 bubbles / kg-day. The choice of CF3I provides optimal sensitivity to both 
spin-dependent (SD) and spin-independent (SI) WIMP couplings, in the first case by its fluorine 
content, in the second via the presence of iodine (I71l2l[73l). This allows a maximally efficient 
exploration of supersymmetric WIMP candidates d?]). 

Bubble production within the chamber is monitored via two trigger mechanisms, the sound 
emission or pressure rise created by their formation and growth, and by the changes that their 
appearance induces in CCD camera images, inspected every ~20 ms. Two orthogonal cameras 
watch the inner volume (Fig. 2), where the active liquid resides within a thin quartz vessel. 
A bellows mechanism balances the pressure difference across the quartz vessel wall ([74]) . The 
precision in the stereoscopic reconstruction of the spatial position of small early bubbles is of 
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the order of their size at the time of the trigger, ~1 mm. A live-time fraction of ~80 % was 
obtained after the adoption of real-time image analysis as the primary trigger. 

The chamber was operated continuously, in a variety of conditions (pressure, temperature, 
presence of radioactive sources, different triggering methods and pressure cycling protocols), 
for a year starting on December 2005, generating an exposure larger than 250 kg-days, with an 
overall mean time between expansions of ~100 s. The emphasis on this first large prototype 
was purely on the functional aspects leading to continuous unattended operation and enhanced 
stability, obtained by implementing the techniques described in ©. In the interest of rapid 
deployment, little attention was paid to alpha-recoil backgrounds such as those induced by 
radon and its progeny (Fig. 1). The data at hand, although entirely dominated by radon-induced 
backgrounds, already yield an improved sensitivity to SD WIMP-proton couplings. 

Calibrations Of importance prior to a WIMP search is an empirical determination of the 
maximum degree of superheat (I75l) achievable before the onset of sensitivity to minimum ion- 
izing particles such as photoelectrons (Fig. 1). This in turn defines the lowest recoil energy that 
can be detected in background-free conditions, and with it the acceptance for a WIMP signal. 
Similarly, it was necessary to demonstrate that sensitivity to low-energy nuclear recoils was 
nevertheless present in those conditions. Two calibration sources were developed for these pur- 
poses. The first one, an intense 13 mCi ^'^''Cs gamma source, was placed next to the chamber's 
outer steel vessel, inside its neutron-moderating polyethylene shield. Fig. 3 shows the response 
of the chamber as a function of operating pressure in the presence of the source and in its ab- 
sence. A Monte Carlo simulation (l6\i was used to calculate the rate of photon-induced electron 
production within the active volume of the chamber. From the difference between this inter- 
action rate and the observed bubble nucleation rate in the presence of the source it is possible 
to obtain a gamma rejection factor (the fraction of interacting gammas inducing bubbles) as a 
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function of degree of superheat. This factor can then be expressed in terms of the calculated 
nuclear recoil energy threshold for each pressure setting (Fig. 3, inset). The minimum-ionizing 
background rejection obtained, ~ 10^^° for a nuclear recoil threshold of ~10 keV, is unmatched 
by any other WIMP detector. For instance, dark matter detection efforts using cryogenic ger- 
manium detectors sensitive to both ionization and heat (Akerib et al. in (1731) ) presently feature a 
rejection factor 10""^ — 10^^ for a similar recoil threshold. As a result of this, moderately super- 
heated bubble chambers are free from any radiopurity constraints in gamma or beta emitters in 
the detector or its neutron moderator shield, including possible elevated rates of ^''C in the active 
liquids. This freedom, when combined with room-temperature operation, translates into large 
gains in construction speed and a reduction in costs. This background insensitivity is intrinsic: 
electron-induced events simply do not take place. 

The second calibration source is a switchable Am/Be neutron source (Fig. 4). It was char- 
acterized using a '^He counter surrounded by moderator. A Monte Carlo simulation (|76l) was 
used to generate the response of the counter, leading to a measured neutron yield of 4.9 n/s, in 
excellent agreement with the (a, n) reaction yield from the ^"^^Am foil sources employed and 
separately characterized. In order to further assess the uncertainty in the very small yield of 
the source (|79|) . the same simulation package and neutron detector were used to characterize 
five commercial neutron sources of known activity. This uncertainty was found to be a modest 
±11%. This switchable source, also placed outside of the steel wall of the recompression cham- 
ber, is moderated by the ~10 cm of recompression fluid around the quartz vessel and is therefore 
expected to produce recoils very similar in energy to those from WIMP interactions (Fig. 4, in- 
set). In comparing the predicted response to this source with actual observations, one must take 
into account additional sources of uncertainty affecting the input to the MCNP-PoliMi Monte 
Carlo (|27]) used to transport neutrons and generate the expected rate and energy distribution 
of recoils in the active liquid. We have appraised these to be dominated by a ±28% from the 
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fiducialization of the chamber geometry and a ±24% from the hardness of the neutron spectrum 
(affected by alpha-particle energy losses on their way to the Be foils (|22)) '). The comparison is 
very good, as can be appreciated in Fig. 4. These predictions were generated prior to inspection 
of the neutron irradiation data. The success of this calibration confirms the predictions for bub- 
ble nucleation thresholds as a function of superheat, sanctioning WIMP limits obtained from 
the device. Switchable sources similar to this one can be used in planned larger chambers for 
periodic studies of response to nuclear recoils, important for instance when looking for WIMP- 
induced modulations in the data (3). To the best of our knowledge this is the first instance of 
their use in a direct search for WIMPs. 

Present limitations of the method Two backgrounds became evident during underground 
operation of the chamber. The first is an excess of bubbles on the wall of the quartz vessel. A 
fraction of their rate may arise from the exposure of the vessel to typical concentrations of radon 
in air before installation. This results in the shallow implantation of its long-lived alpha-emitting 
daughters (|2Q|) . Another source able to explain most to all of the observed surface-alpha rate is 
the ~50 ppb of uranium measured via gamma spectroscopy in the quartz vessel material. An 
origin in alpha emission for these wall events is confirmed by their disappearance at pressures 
higher than 65 (35) psig, the predicted thresholds for alpha-induced bubble nucleation when 
operating at 40 (30) C. Image analysis readily identifies these bubbles as happening on the 
wall, allowing rejection with high confidence. However, since each event is followed by tens of 
seconds of recompression time to ensure recovery of the superheated state in the subsequent de- 
compression (7), an excessive wall rate per unit quartz surface would lead to a considerable dead 
time in larger chambers. A decrease by x 10 ( x 100) in the wall event rate from the presently ob- 
served value is needed to ensure that 0(50) kg chambers under construction sustain a live time 
>60 % (>85 %). In order to achieve this reduction, an etching procedure has been developed 
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in collaboration with the vessel manufacturer, facilitating the removal of implanted daughters 
while preserving surface smoothness and resistance against fracture initiation. In addition, new 
chambers employ synthetic fused silica containing only a few ppt of alpha-emitters. 

The second source of background events from this mechanical prototype detector consists 
of temperature-dependent radon emanations and injections into the chamber, likely from tho- 
riated weld lines and a Viton 0-ring exposed to the inner detector volume. Viton is known to 
release radon at a rate of ~300 mBq/m^ and does not create a strong barrier against its diffusion 
from external sources. Radon decays can lead to single bubbles in the bulk identical to those 
expected from WIMPs (they nevertheless display a very different spectral and temporal depen- 
dence, as discussed below). That the origin of essentially all presently-observed bulk events is 
in radon-associated alphas and alpha-recoils is evidenced by two signatures. First, the distribu- 
tion of times between consecutive bulk events (Fig. 5) displays a clear short-lived component 
corresponding to the the decay sequence — s-^^^Po ^"^^^Pb (both steps are alpha-emitting). 
Two independent analyses of this distribution, including the full decay chain (Fig. 5), the effect 
of recompression (dead) time, fiducial volume cuts and the possibility of other sources being 
present, strongly favor models with ~100% of the events belonging to the ^^^Rn sequence and a 
~100% bubble nucleation efficiency for the combined alpha and alpha-recoil emission. Second, 
the distribution of bulk event rates vs. operating pressure at a fixed temperature is flat below a 
sharply defined onset pressure, the value of which is again in good agreement with theoreti- 
cal expectations (Fig. 5). This is the kind of spectral behavior expected from the response of a 
threshold detector to a mono-energetic source. It has been observed in the past in bubble cham- 
bers intentionally spiked with alpha emitters (|25]) . New larger chambers under construction use 
metallic seals, with all weld lines being non-thoriated, among other precautions against Rn. 
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Improved sensitivity to WIMP couplings 

Even constrained by the background limitations of the mechanical prototype in these early runs, 
an improved sensitivity to the spin-dependent WIMP-proton coupling is obtained from these 
data. This is possible due to the extreme sensitivity to this mode of interaction afforded by the 
large mass fraction of fluorine (29.1%) and the operation of the chamber in conditions where 
only nuclear recoils (alpha, neutron or WIMP-induced) can produce bubbles. Fig. 5 displays the 
characteristic recoil rate expected from example WIMPs, overlaid on the data. Such spectra are 
calculated for different WIMP masses using the theoretical bubble nucleation thresholds (top 
axis in the figures) for fluorine recoils (24). These are in turn used to generate the integrated 
WIMP recoil rate above threshold, following the method and dark matter halo parameters rec- 
ommended in (|25l) . These threshold calculations have been validated by their success in pre- 
dicting the response of the chamber to the Am/Be source (l26l) . Following a procedure similar 
to that in (|27|) . experimental bulk bubble nucleation rates vs. pressure, such as those depicted 
in Fig. 5, are fitted with a model containing two free parameters: the signal from a WIMP of a 
given mass, scalable by a free spin-dependent cross-section, and the response to ^^^Rn alphas, 
expressed as a logistic function with a free overall normalization (Fig. 5). The minimization is 
performed using the MINUIT package (|2^ . its output cross-checked via a standard error ma- 
trix analysis. The best fits favor the null hypothesis alone (response to ^^^Rn alphas only). The 
largest spin-dependent couplings allowed by the data, up to a 90 % C.L., are displayed in the 
form of exclusion plots in Fig. 6. These bounds arise from the weighted average (|25|) of the 
analyses performed on three independent large data sets adding up to 52 kg-days of exposure. 
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Conclusions 



Even prior to any measures against alpha-emitting backgrounds, improved bounds on spin- 
dependent WIMP couplings can be obtained from the application of an old technology, bubble 
chambers, to a new problem, that of dark matter detection. In particular, these new bounds make 
a spin-dependent interpretation pQ| of the DAMA claim of an observed WIMP signal (|37l) very 
difficult. Other experiments have already severely constrained spin-independent explanations. 

In principle, and strictly from the point of view of internal backgrounds, the sensitivity 
of this technique can be improved by as much as six orders of magnitude by reducing alpha- 
emitting contaminants from the uranium and thorium chains to levels similar to the best achieved 
in large neutrino experiments (< 10^^'' g/g) This would guarantee an extensive probing 
of supersymmetric WIMP models along both spin-dependent and -independent couplings (0). 
The possibility exists to add the detection of dark matter to the many contributions already 
made by Bubble Chambers to particle physics. We presently concentrate on the development of 
radon-free chambers with a total target mass of 0(50) kg. 

We gratefully acknowledge the effort and outstanding technical support of the Fermilab 
Staff. This work is supported by NSF CAREER award PHY-02398 1 2, NSF grant PHY-0555472, 
the Indiana University South Bend R&D committee, the Kavli Institute for Cosmological Physics 
through grant NSF PHY-01 14422, and the U.S. Department of Energy. 
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energy (keV) 

Figure 1: Instantaneous stopping power vs. energy for different particles in CF3I, including its 
three recoiling species. Plotted as horizontal and vertical lines are the calculated dE/dx and en- 
ergy thresholds for bubble nucleation at T=40 C and two different operating pressures. Accord- 
ing to the "Hot Spike" nucleation model (8) only radiations in the top right (colored) quadrants 
can lead to bubbles. Notice the absence of this possibility for electrons even toward the end of 
their range, in conditions that nonetheless lead to a sensitivity to recoils of just a few keV, such 
as those expected from WIMP interactions. Alpha particles and their recoiling daughters can 
induce bubble nucleations and the presence of their emitters must therefore be avoided. 
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Figure 2: Families of events in a 1.5 kg CF3I bubble chamber. At high degrees of superheat (~60 
C and atmospheric pressure for this compound), minimum ionizing cosmic ray events reminis- 
cent of those observed in early bubble chambers © are visible (left). At moderate superheats 
(~30 C, 1 atm) the chamber is sensitive strictly to high dE/dx radiation such as nuclear recoils. 
Whereas neutrons can give rise to simultaneous separate bubbles each corresponding to a scat- 
ter (center), WIMPs are expected to produce single bubbles only (right), due to their extremely 
small probability of interaction. The mean free path between scatters is of just a few cm for 
neutrons, leading to an excellent ability to reject them in large chambers. 
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Figure 3: Response of the chamber to an intense Cs gamma source. The expected gamma in- 
teraction rate with the superheated liquid is 3.9x10^ per second, with gamma-induced electron 
energies reaching up to 662 keV. Inset: Intrinsic gamma rejection factor (fraction of interacting 
gammas inducing bubbles) obtained from the exposure to the source (see text). 
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Figure 4: Blind absolute comparison between expected bubble nucleation rate (lines) and ob- 
servations (points) in the active presence of the switchable Am/Be neutron source {top inset). 
The two lines correspond to largely different values of a nucleation parameter "a" (|77l) . the sin- 
gle free factor in the classical theory used to predict bubble nucleation thresholds ®. A fit to 
the data provides an excellent agreement with theoretical expectations (a = 6) (fT^I) . Including 
the uncertainty in the predictions (see text) the same fit simultaneously yields an efficiency in 
the response to this source of 81±^^% (51±i^%) at 40 C (30 C) (errors are 90% confidence 
levels). Calculated thresholds for recoil-induced nucleation are expressed in keV along the top 
axis. Bottom inset: Spectrum of nuclear recoil energies produced by the source (MCNP-PoliMi 
simulation (|27]) ). similar to typical expected WIMP recoil spectra. 
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Figure 5: Top: Distribution of times between consecutive bulk events in the chamber. A model 
(solid line) including the effect of dead time and fiducial volume cuts, based on the triple alpha 
emission from ^^^Rn and progeny, successfully reproduces this distribution (see text). Bottom: 
Distribution of single bulk bubble nucleation rate vs. operating pressure, for two different run- 
ning temperatures. Rates display a flat behavior up to a pressure endpoint. This is characteristic 
of the response to monochromatic alphas and ~100 keV alpha recoils (here from ^^^Rn emana- 
tions, see text). Colored arrows along the bottom axes indicate the predicted onset of sensitivity 
to these particles, in good agreement with observations. As a reference, solid lines correspond 
to the expected signal rate from WIMPs with a mass of 10 and 50 GeV/c^ and a 3 pb cross 
section for their spin-dependent coupling to protons. Also shown is the response function to 
^^^Rn and progeny (dashed lines). Calculated energy thresholds (in keV) for bubble nucleation 
by fluorine recoils are shown along the top axes. 
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Figure 6: Top: Improved limits on spin-dependent (pure) proton-WIMP coupling vs. WIMP 
mass from this experiment (COUPP, the Chicagoland Observatory for Underground Particle 
Physics). Couplings above the line would have produced signals above observed backgrounds 
and are excluded to 90% C.L. Limits from other experiments are also shown (29), as well as the 
(orange) region favored as a possible explanation to an existing claim for WIMP observation 
(130] UTI ). a hypothesis now contradicted by this experiment. Bottom: Similar limits for spin- 
dependent coupling parameters where no assumption is made about the relative strength of the 
coupling to neutrons and protons, but a WIMP mass must be chosen (50 GeV/c^ here) (I32II33I) . 
The region outside of the ellipses is excluded by each experiment. 
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